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Abstract 

Contemporaneous aggregation of individual AR(1) random processes might lead to 
different properties of the limit aggregated time series, in particular, long memory 
(Granger, 1980). We provide a new characterization of the series of autoregressive 
coefficients, which is defined from the Wold representation of the limit of the ag¬ 
gregate stochastic process, in the presence of long-memory features. Especially the 
infinite autoregressive stochastic process dehned by the almost sure representation 
of the aggregate process has a unit root in the presence of the long-memory prop¬ 
erty. Finally we discuss some examples using some well-known probability density 
functions of the autoregressive random parameter in the aggregation literature. 
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1 Introduction 


Aggregation is a critical and widely acknowledged issue in the empirical and theoretical 
literature in economics and other fields. Especially, since the contributions of Granger 
and Morris (1976) and Granger (1980), it is well-known that the contemporaneous ag¬ 
gregation of individual random AR(1) stochastic processes might lead to long memory 
modelsj^ Notably, Granger (1980) considers the case of a Beta-distribution for the random 
autoregressive parameter and thus points out that the long memory property depends 
on the behavior of the density of the random autoregressive parameter near unity and 
that the common and idiosyncratic components might exhibit a different degree of long 
memory. In particular, Zaffaroni (2004) generalizes these results by studying the limit 
of the aggregate process with a quite flexible (semi-parametric) assumption regarding 
the behavior near unity of the probability density function of the random autoregressive 
coefficient and makes clear the asymptotics of both the common and idiosyncratic parts 
of the aggregated process. Among others results, Zaffaroni (2004) shows formally that 
the more concentrated is the distribution of the random autoregressive coefficient near 
the unit, the stronger is the long-memory property of the limit aggregated processR Fol¬ 
lowing these contributions, we study the aggregation of heterogenous individual AR(1) 
stochastic processes that leads to the long-memory property. In contrast to the litera¬ 
ture, we focus on the infinite autoregressive representation of the limit of the aggregate 
process, and especially the sum of the autoregressive coefficients, rather than the usual 
infinite moving representation. 

Indeed the use of the (equivalent) infinite autoregressive representation of the limit of the 
aggregated process might be meaningful in various contexts of the aggregation (panel) 
literature. As pointed out by Lewbel (1994), the autoregressive representation is insight¬ 
ful among others when estimating the (macro) aggregate dynamics with unobservable 
individual series in the presence of unobserved parameter heterogeneity, when identifying 
and estimating certain distributional features of the micro parameters from aggregate re- 

^In the sequel, a stochastic process is said to have a long memory property if its autocovariance 
function is not summable (Beran, 1994; Beran et ah, 2013). 

^For a recent survey, see Leipus et al. (2013). 


2 



lations (disaggregation problem), or when deriving the aggregate dynamics with nnknown 
common error termsj^ Obviously, if the limit of the aggregate process is represented al¬ 
most surely by a short memory process, the inhnite autoregressive representation and 
especially the determination of its persistence (e.g., through the sum of the autoregres¬ 
sive coefficients) easily obtains using standard results of time series analysis (Brockwell 
and Davis, 2002). In contrast, if the limit of the aggregate process is represented almost 
surely by a long memory process, the convergence and thus the determination of the 
series of autoregressive coefficients is challenging since the series of moving average co¬ 
efficients is no longer absolutely summable. Such a characterization might be necessary 
for instance in the case of the estimation of the aggregate dynamics with unobservable 
individual series in which hnite parameter approximation for the inhnite lag distribution 
is required—the autoregressive representation of the limit aggregated process displays an 
inhnite parametrization whereas the econometrician has only hnitely many observations 
and thus hnite-parameter approximations might lead to near observational aggregate 
processes with diherent statistical properties as long as the sum of the autoregressive 
coefficients is not correctly identihed. In this respect, our paper tackles this issue and 
proposes a new characterization of the autoregressive representation of limit long-memory 
aggregate processes that results from individual random AR(1) stochastic processes. 

Notably aggregation of individual random AR(1) processes is analyzed in this paper under 
the assumptions that some common factors exist, the (positive) random autoregressive 
coefficient takes values in [0,1) with probability distribution /i, and all of the noncentral 
moments exist. In particular, the series of all of the noncentral moments might be either 
absolutely summable or divergent and the limit of the aggregated stochastic process sat- 
ishes almost surely both an inhnite moving average and autoregressive representations. 
Within this framework, we show that the sum of the autoregressive coefficients of the 
stochastic process dehned from the limit in C? of the aggregate process, say ^ a^, might 

fc>i 

equal one (respectively, less than one) when the limit aggregate process has the long 
memory property (respectively, short memory property). Say differently, the divergence 
of the series of noncentral moments is fully equivalent to the presence of a unit root in 
^See also Pesaran (2003), Pesaran and Chudik (2014), Jondeau and Pelgrin (2014a,b). 
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the representation of the stochastic process dehned by the limit of the aggregate process. 
In so doing, we consider the mapping between the moving average coefficients and the 
autoregressive coefficients of the limit aggregate process and make use of complex analysis 
when the series of the moving average coefficients is divergent. Complex analysis is called 
for at least two arguments. On the one hand, the existence of the limit linir^i- a{r) of the 
(ordinary) generating function of the sequence of the autoregressive coefficients does not 
insure its equality to even the convergence of this sum. On the other hand, we 

cannot apply a standard Hardy-Littlewood Tauberian theorem since, to the best of our 
knowledge, there is no general proof of the positiveness of the autoregressive coefficients 
for all probability distributions jj, dehned on [0,1). Interestingly such a complex analy¬ 
sis might be used to study the behavior on the unit circle of an inhnite parameterized 
stochastic process, which is dehned (almost surely) in but does not belong to . 

The rest of the paper is organized as follows. In Section 2, we discuss the main assump¬ 
tions regarding the individual random AR(1) stochastic processes and then we derive 
the limit aggregate process. In Section 3, we provide the two main results of our paper, 
namely the determination of the sum of the autoregressive parameters and the charac¬ 
terization of the stochastic process dehned by the limit in of the aggregated process. 
In Section 4, we assume that the distribution of the random autoregressive coefficient is 
subsequently a Beta distribution (of the hrst kind), a uniform distribution and a random 
polynomial density function. In so doing, we apply the results of Section 3 and generalizes 
some results in the literature. Proofs are gathered in the Appendix. 

2 Aggregating individual AR(1) processes 

In this section, we hrst discuss the assumptions regarding the individual AR(1) processes. 
Then we derive the limit aggregate process. 
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2.1 Individual dynamics 

Consider the individual random AR(1) model for i = 1, • • • , N and t G zQ 

"h (f) 

where (pi denotes an individual-specihc (random) parameter drawn from a hxed random 
variable p and Vi^t is an error term that can decomposed into a common component, e^, 
and an idiosyncratic (individual-specihc) component, 

Vi,t = + Tji^f ( 2 ) 


The macro variable results from the aggregation of micro-units, with the use of time- 


N 


invariant nonrandom weights Wn = (tci, • • • ,^c^^)^ with '^Wi = 1, so that the aggregate 


2=1 


N 


variable is dehned as the weighted average of the micro-units The fol- 

2=1 

lowing assumptions hold: 


Assumption 1: pa. hxed random variable with probability distribution /i the sup¬ 
port of which is in [0,1]. 

Assumption 2: The moments Uk = E(y9^) exist for all integer k>l. 

Assumption 3: e* and are white noise processes with means of zero and variance of 
CTg and cr^, respectively; et and 7]i^t are mutually orthogonal at any lag and lead. 
Assumption 4: Realizations of p are independent of those of e^. 

Assumption 5: As A ^ cx), lllTjvll = O (A-V2) and Wi/ ||hRv|| = O for all 

i e N. 


Before discussing our main assumptions, we introduce some notations. Let T-Lx de¬ 
note the Hilbert space generated by all of the random variables that compose the panel 

^Alternatively one might assume that t S N and consider an asymptotically stationary limit aggregate 
process (Gongalves and Gourieroux, 1994). 

®Such dynamics have been used in economics (among others) to represent consumption expenditures 
across households (Lewbel, 1994), consumer price inflation across subindices (Altissimo et ah, 2009), real 
exchange rates across sectors (Imbs et ah, 2005), or real marginal cost across industries (Imbs et ah, 
2011 ). 
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'Hxi the Hilbert space generated by the stochastic process 'Hx,t 

and 'Hxi,t the corresponding snbspaces of Hx and Hxi up to time t. Assnmption 1 ontlines 


0 


that is a random variable with distribntion on the snpport [0,1) snch that (p G f](BT-Lxi,t 

t * 

This assnmption is consistent with many parametric specihcations of the cross-sectional 
distribntion of p. We only rnle ont sitnations in which some individnal processes are not 
almost snrely (asymptotically) stationary, i.e. P (|<p| > 1) > 0 (see farther). In particular, 
Assumption 1 includes the Beta distribution, B{p,q). In this case, the representation of 
the (limit) aggregate process does depend on the mass distribution of the Beta distri¬ 
bution around unity: the smaller is q, the larger is the mass of the distribution around 
unity (Zaffaroni, 2004). In contrast, imposing the condition 0 < p < c < 1 for some 
constant c will guarantee that there are no individual unit root parameters that would 
dominate at the aggregate level (Zaffaroni, 2004) and that the limit aggregate process (as 
N —)■ cxo) displays short memory with an exponentially decaying autocorrelation func¬ 
tion. Assumption 2 insures that noncentral moments Uk = E(v9^) of any (nondegenerate) 
random variable p, dehned on [0,1), satisfy: 1 > > • • • > ma: > 0, V/c > 1, and ma, — )■ 0 

as /c —)■ cxo. 


Assumption 3 dehnes the statistical properties of the two components of the error term, 
the common shock (factor) and the idiosyncratic shock. Several points are worth com¬ 
menting. First, without loss of generality, one might assume that the stochastic process 
,N,tez is weakly linearly exchangeable (Aldous, 1980), i.e. is 

a purely non deterministic (regular) covariance stationary process and the covariance 
structure is invariant by translation with respect to the time index and invariant by 
permutation with respect to the unit index. Second, the stochastic process {et)tez is as¬ 
sumed to be known and given Given that T-Lx = (BTixi and thus Tix^t = ®'Hxi,t, one has 

i i 

et G 'Hx,tf]{'Hx,t-i)±- Third, taking that Tixi = B-e (B'Hxi-e, the idiosyncratic stochastic 
process is such that G 'Hxi-e,tr\i'^xi-E,t-i)±- Fourth, the assumption that the vari- 


®Note that the set (—1,0) has been excluded from the support of y for simplicity’s sake. 

^Without altering our main results but at the expense of further assumptions, (i) the common error 
term at time t might be multiplied by the realization of a (scaling) random variable k to introduce some 
form of heteroscedasticity, (ii) the stochastic process (et)tgz might be assumed to be unknown (using 
a two-step procedure), and (iii) multiple independent common error terms might be introduced in the 
specification (Zaffaroni, 2004). 
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ance of the idiosyncratic shocks is the same across individuals might be too stringent in 
practical examples (e.g., in economics or in hnance) and can be generalized by assuming 
heteroscedastic idiosyncratic error terms. Eq. (|2]) together with Assumption 3 provide 
a parsimonious form of (statistical) cross-sectional dependence, which is common in the 
aggregation literature (Forni and Lippi, 2001; Zaffaroni, 2004). Obviously one may in¬ 
troduce dependence on the idiosyncratic part through cross-correlation and/or spatial 
correlation (Chudik, Pesaran and Tosetti, 2011; Pesaran and Tosetti, 2009). Such corre¬ 
lation has no effect on our results as long as it is sufficiently weak (e.g., local correlation) 
so that the error structure belongs to the class of approximate factor structure models 
(Chamberlain and Rothschild, 1983) jj. Assumption 4 avoids any correlation between the 
error terms and ip. Assumption 5 is a granularity condition, which insures that the weights 
used to dehne the aggregate process are not dominated by a few of the cross-sectional 
units (Gabaix, 2011; Pesaran and Chudik, 2014)]^ 


2.2 Aggregate Dynamics 

The empirical cross-sectional moments of ip are Ejv = '^Wiip^, \/k > 1. For sake of 

i=l 

simplicity and without loss of generalization with respect to Assumption 4, we assume 
that Wi = 1/N for all i. Consequently, as A^ —)■ cxo, Ea? ( 9 ?^) Uk- 


Using Eqs. (HD-da), the exact aggregate dynamics can be written asl/j 


N ^ V 

n (1 - p>jL) XN,t = n 

j=l i=l j^i 

or equivalently 


(3) 


+ ^^(1-<^iL) Sm (4) 

i=l i=l 

®The condition is that the maximum eigenvalue of the covariance matrix of % = (?7i,t,-'- jijN.tY 
remains bounded as the cross-section dimension increases. 

®Our results extend to the case of (time-varying) stochastic weights. Such an extension requires at 
least that the weights be distributed independently from the stochastic process defining the random 
variable. 

^°Put differently, it is an ARMA(N,N-1) in the absence of common roots in the individual processes 
(Granger and Morris, 1976). 
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where L is the lag operator {zt-i = Lzt). Taking Eq. (jl]), we can characterize the 
asymptotic behavior of both the idiosyncratic component and the common component. 
This is done in the following proposition. Results are known but are reported here for 
sake of completeness. 

Proposition 1 Suppose that Assumptions 1-5 hold. Given the disaggregate model de¬ 
fined in Egs. m-0i, the limit in C? of the aggregated process as N ^ oo satisfies (almost 
surely) the two eguivalent representations: 

OO 

Xt = '^UkCt-k (MA form), (5) 

k=0 

OO 

Xt = akXt-k + et (AR form), (6) 

k=l 

and Eat = E as N ^ OO. The seguence {ak,k > 1} 

satisfies the recurrence relation : 

Ai = g) , Ak+i = {Ak - Ok) (f (7) 

Proof: Gongalves and Gourieroux (1988), Lewbel (1994). 

Several points are worth commenting. First, as shown by Robinson (1978), with the 
exception of a degenerate distribution for (p (e.g., Dirac distribution), the dynamics of 
the limit aggregate process is richer than the individual dynamics because of the non- 
ergodicity of the individual random AR(1) process. Second, using the inhnite moving 
average representation and the positiveness of the moments, the (limit) aggregate pro¬ 
cess displays short memory if < +oo whereas it has a long-memory property 

if Er= ,oUk = -l-cxo. It is worth noting that (i) the (limit) aggregate process is in 
but not necessarily in C^, and (ii) the sum of the autoregressive coefficients is (abso¬ 
lutely) convergent if the aggregate process has short memory (since the spectral repre¬ 
sentation is unique). Third, the central moments of the cross-sectional distribution of 
ip can be easily obtained from the inhnite autoregressive representation of the aggregate 
process. This is useful when considering the standard disaggregation problem in statis- 


where X^g —t Xt 
where = E [Ak] 


tics. For instance, the first four cross-sectional moments are E [<yc] = oi, V [<yc] = 02 , 
S [if] = (03 — 0102) / (02)^^^, and K [ip] = (04 — 20103 -|- 0^02 -|- a ^) / (02)^- Fourth, Equa¬ 
tion ([n]) shows that aggregation leads to an inhnite autoregressive model for Xf (see 
Robinson, 1978, Lewbel, 1994). Notably, using Eq. ([7]), the autoregressive parameters o^ 
are nonlinear transformations of the noncentral moments of ip and satisfy the following 
non-homogenous difference equations (for k > 2): 

k 

ttk+l — E [ylfc_|_i] ^ ^ OirUk—r+l 

r=l 

ai = E [p]. 


Fifth, the persistence of the stochastic process dehned by the limit of the aggregated 
process can be dehned as the sum of the autoregressive coefficients, which is denoted 
a(l) = ^k- More specihcally, if the limit of the aggregate process belongs to and 

thus (Xt)iez is a short-memory stochastic process, then it is straightforward to show that 


a(l) 



-1 


E 


l-p 


or equivalently 


( 1 -«(!))■' 



In particular, the limit of the aggregate stochastic process has short memory if and only 


if E 


l-^p 


< 00 (Robinson, 1978; Gongalves and Gourieroux, 1994). In the spirit of 


Zaffaroni (2004), a sufficient condition on the probability density function h oi p for 
the short memory property is that there exists a G (0,1) and a constant C such that 
lim ,, = C. Sixth, we can dehne the generating function of the terms as follows. 

3.^1- 

Taking the recurrence relation, one can write formally 




zp 

1 — zp 
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and thus 


akz’^ = UkZ^ - ttkZ^ Ukz'". 
k>l k>l k>l k>l 

Let m denote 

m{z) = '^UkZ^, 

k>l 


we obtain the formal generating function of the terms 


= Y1 

k>l 


ttkZ^ 


m{z) 

1 + m{z) ’ 


( 8 ) 


and thus the mapping between the generating function of the inhnite autoregressive lag 

polynomial and the one of the inhnite moving average representation. As explained in 

Section 3, Eq. (jH]) is fundamental since our proof is intimately related to the function 
m 

1 + m 


3 Aggregate long-memory process and the sum of 
the autoregressive parameters 

In this section, we show that the sum of the autoregressive coefficients equals one in the 
presence of a (limit) aggregate process in but not in In so doing, we emphasize that 
complex analysis is required to obtain the convergence of the series Ylik>i^k- Then we 
provide a characterization of the (limit) aggregate long-memory process in the presence 
of individual random AR(1) processes. 


3.1 The function m 

Before showing the convergence of the series of autoregressive parameters, we need an 
intermediate result regarding the m function, and especially to make clear the relationship 
between the sum of the moving average coefficients '^k>i'^k and m{r) when r —)• 1“. 
Proposition [2] clarifies this link and turns to be extremely useful when characterizing 
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the infinite sum of the autoregressive coefficients of a (limit) aggregate long-memory 
process m 

OO 

Proposition 2 = +oo if and only if lim ,,^. 1 - m(r) = -|-cx) 

k=l 

Proof : See Appendix A. 


3.2 The series a(z) = 


m{z) 


and the convergence of Y2k>i 


1 + m{z) 

One key issue to study the convergence of the autoregressive coefficients is that the ex¬ 
istence of the limit, lim^^i- a(r), does not insure its equality to even the 

convergence of this sumo Supplementary Tauberian conditions are needed for this kind 
of results (Hardy, 1949; Titchmarsh, 1939; Korevaar, 2004). In particular, the conver¬ 
gence of this series has to be studied by making use of complex analysis, especially in the 
case of a (limit) aggregate long-memory process. 


All of the proofs related to Theorem 1 are gathered in Appendix B. To summarize, 
the proof of the convergence of the autoregressive coefficients proceeds as follows. We 
first define m as an analytic function of the complex variable z within the open disc 
around 0 of radius one, 11(0,1), and rewrite m{z) as an integral on [0,1) with parameter 


z G Z1(0,1). This allows showing that it can be continuously extended to Z1(0,1)\{1} 
(see Lemma 1 in Appendix B). Then a second lemma (see Lemma 2 in Appendix B) 
proves that the function a{z) = is well defined in the disc /1(0,1) (i.e., the de¬ 

nominator does not vanish) and then it can be extended to the whole closed disc /1(0,1). 
For this purpose, two cases must be studied, according to the nature of the series 
if ^ Mfc < C )0 the function m is continuous in the closed disc, so is a ; and if ^ Ufc = oo 
then \m{z)\ —?■ oo when z ^ 1 in 11(0,1), therefore a{z) 1 when z —)■ 1 in Z1(0,1). 
Finally, a third lemma shows that the analytic function a is univalent and provides the 


^^Note obviously that = -too if lim ,—>. 1 - ra(r) = -too. 

k=0 

OO 

^^For instance, consider the power series defined for |r| < 1, h{r) = ^ (—Then h{r) = and 

k^O 

oo 

lim h{r) = i. However X](“l)^ convergent. 

r-s-l- 0 
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use of a Tauberian condition (see Lemma 3 in Appendix B). This then allows to prove 
the convergence of the series '^an (Theorem 1). 


Theorem 1 Let {a^, k > 1} denote the sequence defined in Proposition [H 

OO 

is convergent and 

k=l 


+00 , N 

m(r) 

> Ofc = hm - 

1 + m{r) 

k=l ^ ’ 


The series 


Proof : See Appendix B. 


Taking Proposition [2] and Theorem [H it is then straightforward to determine the sum of 
the autoregressive coefficients. 


Proposition 3 Let {ak,k > 1} denote the sequence defined in PropositionUi The sum 

OO 

of the autoregressive coefficients, equals one if and only z/lim^^.!-m(r) = +cxo or 

k=l 

OO 

equivalently if and only if = +cxo. 

k=l 

Three points are worth commenting. First, if the limit aggregate process is a second- 
order stationary process (Proposition [T]) and the series of its moving average coefficients 
is absolutely summable, then the sum of the autoregressive coefficients is less than one. 
Notably, this result obtains with classical time series results whereas there is a need of 
complex analysis when the series of the moving average coefficients is not absolutely 
summable and one studies what happens at the pole z = 1. A second and related point is 
that the behavior of the series of the autoregressive coefficients depends on whether the 
limit second-order stationary aggregate process belongs to or not. Consequently, as 
explained below, this provides a new characterization of a (limit) aggregate long-memory 
process as a result of the aggregation of random AR(1) processes. Third, as stated in 


Corollary [T], the function $( 2 ;) = 1 — a{z) admits only one zero on £>(0,1). 


Corollary 1 Let $( 2 ;) = 1 — a{z). Then z = 1 is a zero of the function <F, which is 

_ OO 

defined on D{0, 1), if and only if = + 00 . 

k=l 
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Proof : This is a straightforward implication of Theorem [T] 


Corollary [T] establishes that the representation of the stochastic process dehned by the 
limit of the aggregate process admits a unit root whereas the aggregate pro¬ 

cess is weakly stationary. There is a one-to-one relationship between the long memory 
property and the presence of a unit root of the inhnite autoregressive limit aggregate 
process. 


3.3 Time series implications 

Taking Theorem 1 and Proposition [3l we are now in a position to provide a new charac¬ 
terization of a (limit) aggregate long-memory process that results from the aggregation 
of individual random AR(1) processes. 


Theorem 2 Suppose that Assumptions 1 and 2 hold true. Let {Xt,t G Z) denote a 
long-memory process with the following Wold decomposition 

OO 

k=0 


OO 

where the Uk terms are nonnegative and ^Uk = oo. Then the aj terms of the equivalent 

k=0 

infinite autoregressive representation, 

OO 

Xt = E Q-kXt-k + Q, 

k=l 


which are defined from a^+i 


^fc+i r '.,—1 r+i and tti 


Ui, satisfy flfc 


1 . 


Several points are worth discussing. First, the (limit) aggregate long-memory process 
does not belong to the class of ARFIMA or self-similar stochastic processeso In partic¬ 
ular, it is not possible to express the moving average coefficients of a fractional ARIMA 
process (by binomial expansion) such that they match the moving average weights of the 


continuous-time stochastic process (R) is said to be self-similar with self-similarity paramter H, if 
for any sequence of time points ti, - ■ ■ ,tk and any positive constant a, c~^{Yati , • ■ • , Wtk) the same 
distribution as (R^, • • • , R^). 
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long-memory process defined in Theorem [2l In the same respect, matching the aggregate 
moving average coefficients of the (limit) aggregate process reqnires non-constant values 
of the self-similarity parameter H (Beran, 1994; Beran et ah, 2013). Second, the fact 
that the sum of the autoregressive coefficients of the (limit) aggregate process is equal 
to one is also consistent with the standard dehnition of long-memory processes (Beran, 
1994; Beran et ah, 2013), i.e. there exists a real number a G (0,1) and a constant Cp 
such that lim = 1 where p{k) = 7x(^)/7x(0) (with 7 x(h) = E /{I — ’~p^)\) is the 
autocorrelation of order k of the (limit) aggregate stochastic process (Xt). 


Third, as pointed out by Beran (1994), observing long-range dependence in an aggre¬ 
gate time series (e.g., at the macro-level) does not necessarily mean that this is due 
to the genuine occurrence of long memory in the individual series (micro-level). This 
might be induced artihcially by aggregation. Say differently, identifying the source of 
long memory would require to look carefully at the behavior of the possibly unobservable 
individual series. Fourth, Theorem |2] has some implications, which are beyond the scope 
of this paper, regarding some aggregation problems often encountered in the theoretical 
and empirical aggregation research (Pesaran and Chudik, 2014), namely the derivation 
of the macro dynamics from heterogenous individual dynamics and the identihcation and 
estimation of certain distributional features of the micro parameters from aggregate rela¬ 
tions (disaggregation problem) o For instance, in the case of the disaggregation problem 
when individual series are not available, since the autoregressive specihcation in Proposi¬ 
tion [T] displays an inhnite parametrization and the econometrician has only hnitely many 
observations, one might proceed with a hnite-parameter approximation for the inhnite 
lag distribution (Sims, 1971, 1972; Faust and Pepper, 1997; Potscher, 2002) and might 
account possibly for the constraint a(l) = 1. Notably, hnite-parameter approximations 
might lead to near observational aggregate processes with diherent statistical properties 
and thus to the incorrect identihcation of distributional features of the micro parame- 
terso Other applications might concern, among others, the estimation of the aggregate 


Lewbel (1994), and Jondeau and Pelgrin ( 2014b). 

^^The convergence of the estimates of the approximation is not sufficient to guarantee the conver¬ 
gence of some functions of those parameters—pointwise convergence does not imply (locally) uniform 
convergence. 


14 








dynamics in the presence of unobserved heterogeneity when individual series are not 
available, the derivation of the macrodynamics in heterogeneous agents models, or the 
reconciliation of micro (panel) evidence and macro facts (e.g., the purchasing power par¬ 
ity, the nominal price rigidity). 


Finally it is worth noting that Theorem |2] applies in a broader context than the ag¬ 
gregation of random AR(1) stochastic processes. Indeed, any stochastic process dehned 
(almost surely) in C? that admits a Wold decomposition with decreasing and nonnega¬ 
tive moment moving average coefficients and does not belong to displays a unit root. 
This comes from the formal identity between the generating functions of (a^) and {uk), 

'^k>l — '^k>l '^k^^ ■ 

4 Examples 

In this section, we review three examples, namely the Beta B {p, q) distribution (of the 
hrst kind), the uniform distribution with p = q = 1, and the random polynomial density. 


Beta distribution Following Morris and Granger (1976), Granger (1980) and Gongalves 
and Gourieroux (1994), we assume that p is Beta-distributedD£ 


R(p, g; x) = ^l[o,i)(a::), P > 0, g > 0. 

In this respect, to the best of our knowledge, PropositionSlprovides a new characterization 
of the series of autoregressive coefficients (i.e., the persistence of the (limit) aggregate 
process (W))- 

Proposition 4 Suppose that p is Beta-distributed (of the first kind). Given the disag¬ 
gregate model defined in Eqs. (QP-ID, the series of the autoregressive coefficients of the 
limit aggregate process defined in Proposition [I] is given by: 

OO 

- Ifq > 1, then 

k=l 

^®Note that Gongalves and Gourieroux (1994) study extensively the aggregation of individual AR(1) 
processes in which ip is Beta-distributed (after an homothety) and provide a discussion regarding aggre¬ 
gate long-memory processes. However they do not consider the sum of the autoregressive coefficients. 
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- Ifq<l, then J2^k = 


k=l 


Proposition 0] can be shown as follows. Taking the integral form of m{r) we have 

1 + m r = / -^- dx, 

r(p)r(g) Jo 1-rx 

and thus (by the monotone convergence theorem) 

lim (1 + m(r)) = [ x^-^il - x^-^dx. 

>'^1- r(p)r(g) Jo 

Consequently, this is a convergent integral if and only if g > 1. In this case, 


r n ^ r(p + g) r(p)r(g-1) p 
r(p)r(g) r(p + g -1) g-l 


and 


+ 00 


E ttk = lim 

r^^ - 


m r 


p 


k=l 


r-i-1- 1 + m{r) p + q — 1 


On the other hand, it follows that lim + m{r)) = +cxo if and only if g < 1, and 

ZlS = 1 (Theorem 1). 


Uniform distribution We now assume that p = g = 1 such that the random variable 
(p is uniformly distributed over the interval [0,1) and non central moments are given by 

^ FTT- 

Proposition 5 Suppose that p is uniformly distributed over the interval [0;1). Given 
the disaggregate model defined in Eqs. m-m, the autoregressive coefficients of the limit 
aggregate process defined in Proposition [I] are given by: 



where Ik = x(x — l) ■ ■ ■ (x — k + l)dx has the same sign as (—1)^ ^ for k > 1. Moreover, 

E OO 1 

k=i — 1 - 

Proposition [5] can be shown by using either Theorem 1 or a new lemma provided in Ap¬ 
pendix C. Notably, the coefficients of the series a{z) = m(^)/(l-|-m(.s)) can be computed 
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as follows. First, the generating moment series is 


1 + m{z) 




E 

fc >0 


k + 1 


log(l - z) 
z 


where < 1 so that 


/ \ miz) 
a{z) = , . . . = 1 + 


1 + m(z) log(l — z) 

Second the expression of is derived by a power series development of the function 

a. Indeed, using the notation 'ip^z) = log(l — z), one has 


f ^ ^-- 

/o log(l - Z) 


where the power series development of = (1 — zY is dehned to be: 


(1 - 2)^ = ^ P J (-1)"^" = 1 + 5^ ^-^z^x{x - 1) • • • (x - n + 1). 


n>0 


n>l 


Since this series converges absolutely for \z\ < 1 and uniformly for x in [0,1], one obtains 


'1 

x-ipiz)^ _ 


^dx = / (1 — zYdx = 1 + ^2 


-1)^ 


■z^L 


n>l 


n\ 


where /„ = x(x — 1) • • • (x — n + l)dx for n > 1. Finally, 


/ e^^^^^dx = —^-- 

fo log(l - 




k>l 


and thus 

= E 

k>l 

Several points are worth commenting. On the one hand, the non-negativeness of the 
autoregressive coefficients for the uniform distribution allows for the use of a standard 
Hardy-Littlewood Tauberian result and thus might not require Theorem 1. However 
the non-negativeness of the autoregressive coefficients is not proved for all probability 
distributions p with support [0,1). On the other hand, using Theorem 1, since the 
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moving average coefficients of the (limit) aggregate process constitute a harmonic series, 
the corresponding series diverges and thus the sum of the autoregressive coefficients (of the 
limit aggregate process) equals one. Finally, Proposition [S] extends the result of Linden 
(1999) in which the behavior of the aggregate process is studied with the autocorrelation 
function. 


Random polynomial density We consider as a last example the case of the poly¬ 
nomial aggregated AR(1) model. More specihcally, we suppose that has a continuous 
distribution over [0,1] that can be represented by a polynomial of degree d > 1 : 

d 

s=0 

where ~ ^ density has to be integrated to one) and / is non-negative in 

[0,1]. One key feature of the polynomial density function relative to the Beta distribution 
is that it can be multi-modal (for polynomial density of third order or above). Similarly 
to the Beta distribution, two cases are considered in Proposition [HI the sum of the 
autoregressive coefficients does depend on whether /(I) = 0 or /(I) > 0 with /(I) = 
Co + Cl -|- • • • -f Crf is the value of the density at x = 1. 

Proposition 6 Suppose that (p has a polynomial density function of order d. Given the 
disaggregate model defined in Egs. m-m, the series of the autoregressive coefficients of 
the limit aggregate process defined in Proposition [I] is given by: 

- ///(I) = 0, then 

OO 

a(l) = 

k=l 


= 1 - 


d—1 n 

EE 

n=0k=0 


^n — k 
n+1 


CX) 

- If f{l) > 0, then 


k=l 


1 . 
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Indeed, starting from the polynomial density function of order d, the non-central moments 
Uk = E k>0, are given by: 


E [P] = 


s=0 


s + k + 


and thus the generating moment series, 


1 -|- m{z) 


k>0 


oo d 

y(Y —— 

^ s -|- A: -|- 1 

k=0 s=0 




is convergent at least for \z\ < 1. Therefore, one needs to study the convergence of the 
series 


1 + m(l) = 


fc =0 s =0 


s + k + 1 


K d 

lim 

K^oo 

k=0 s=0 



Cs 

s + k + 1 


The terms of the double sum S{K) = J2k=o array with K + 1 rows 

and d -|- 1 columns. Without loss of generality, we can suppose K > d. Let n = s -f- fc, we 
have n — d < k < n. Taking diagonal sums along the lines s + k = 0,1,... ,d + K of the 
array, one can write S{K) = Si{K) + S 2 {K) + S 3 {K) with 


S,{K) 

S2{K) 


d n 

EE 

n=0 k=0 
K 


E E 


Cn-k , Co -|- Cl Cq + Cl + ■ ■ ■ + Cd 

— = <=» + ^- + '" +-JVI- 

K ^ n K 

E ETT E = /(I) E 


n -I- 1 

n 

^n—k 


n + 1 n + 1 

n=d-\-l k=n—d n=d-\-l k=n—d 


n=d-\-l 


n + 1 


and 


d-\-K n d-\-K ^ n 

ett ^ 

n=K-\-l k=n—d n=K-\-l k=n—d 

_ Cl + C2 + ■ ■ ■ + Cd C2 + ■ ■ ■ + Cd _ Cd 

K + 2 K + 3 K + D + l 


The hrst sum, Si{K), is hnite and independent of K. The second sum, S 2 {K), is clearly 
convergent as iL —?■ oo if and only if /(I) = 0. Finally, the third sum, S 3 {K), is hnite 
and its limit is zero when K approaches inhnity. In this respect, the generating moment 
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series 1 + m(l) converges if and only if /(I) = 0. In this case its limit is Si = Si{K) with 


Si{K) 


, Co + Cl ^ , 

Co H-z-[-■■■ + 


Co + Cl + ■ • • + Crf 

d + 1 


In contrast, when /(I) > 0 the series J2k>o'^k diverges and thus the sum of the au¬ 
toregressive coefficients equals one. Consequently, the (limit) aggregate process displays 
long-range dependance. 


In this respect. Proposition [6] provides formally the sufficient condition discussed by 
Chong (2006), i.e. /(I) > 0 is sufficient to establish the long memory properties of the 
(limit) aggregate process. Moreover, if we assume that the polynomial density function 
is of order zero and cq = 1 , we end up with the uniform distribution (with /(I) > 0) 
and thus Proposition [5l On the other hand, if we assume that the polynomial density 
function is of the form ax{l — x) for a 7 ^ 0 and x G [0,1), then /(I) = 0 and the long-run 
persistence is given by a(l). In contrast, if the probability density function has a support 
on [0;1], the aggregation process leads to a generalized integrated process (Lin, 1991; 
Granger and Ding, 1996). As a hnal remark, it is worth emphasizing that any distribu¬ 
tion such that the generating function of the sequence of the moments is not convergent 
leads to a long-memory process characterized by Theorem 2. 


5 Conclusion 

In this paper, we study the aggregation of individual random AR(1) processes under 
the assumptions that some common factors exist, the (positive) random autoregressive 
coefficient takes values in [0,1) with probability distribution p, and all of the noncen¬ 
tral moments exist. Notably we show by making use of complex analysis that sum of 
the autoregressive coefficients equals one in the presence of limit aggregate long mem¬ 
ory processes: the divergence of the series of noncentral moments is fully equivalent to 
the presence of a unit root in the autoregressive representation of the stochastic process 
dehned by the limit of the aggregate process. We then illustrate our results using some 
prominent examples of distribution for the aggregation of random autoregressive AR(1) 
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processes. This provides some new insights that might deserve some empirical applica¬ 
tions and some theoretical developments, as for instance the disaggregation problem. 


21 


Appendix: Proofs 


Appendix A 

Proof of Proposition [2J Taking that the sequence of the moments 

Uk = = [ x^dli{x) 

is positive and decreasing, one has hmfc_j.+ooMA; = 0 (by monotone convergence applied 
to the sequence {x^)k>i and thus the radius of convergence of the series 
least 1. 

Let (r„) an increasing sequence in [0,1) such that hm„_j.+oo = 1- For x E [0,1) the 

T X 

functions fnix) = — - -are positive therefore by Fatou’s Lemma we get 

1 - r^x 

CP C “l“00 

lim inf m(r„) = lim inf / /ndp(x) > / liminf /„(ip,(a;) = / - d^{x) = y Uk 

” ” ^[ 0 , 1 ) J[o,i) " i[o,i) 1-2; ^ 

Thus if '^Uk is divergent then m{r) —)■ +cxo when r —)■ 1“. Moreover by the Abel 
Theorem (Titchmarsh, 1939, pp.9-10) if '^Uk is convergent then linv_j.i- m(r) is hnite 
and linv^i- m(r) = 'Y^Uk- The equivalence then follows. 

Appendix B 

In Appendix B, we provide the proof of Theorem 1. In so doing, we proceed with three 
lemmas. Notably we hrst dehne m as an analytic function of the complex variable z in 
the disc /1(0,1) and rewrite m{z) as an integral on [0,1) with parameter G hl(0,1) 
(Lemma 1). We then show that it can be continuously extended to /1(0,1)\{1}. Then 
Lemma 2 proves that the function a{z) = Ym{z) dehned in the disc D{0, 1) (i.e., 

the denominator does not vanish) so that it can be extended to the whole closed disc 
11(0,1). In this respect, two cases must be studied, according to the nature of the series 
Y^k- ii Y'^k < oo the function m is continuous in the closed disc, so is a ; and if 
Y'^k = oo then \m{z)\ —oo when z —)■ 1 in 11(0,1), therefore a{z) 1 when z ^ 1 in 
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D(0,1). Finally, Lemma 3 proves that the analytic function a is shown to be univalent 
and provides the use of a Tauberian condition. Then Theorem 1 is proven. 

Let m the function defined in the open disc £*(0,1) = {^ G C with | 2 r| < 1} by 

+ 00 

m{z) = 

n=\ 


This function is analytic. Moreover, by a classical theorem, on the boundary (7(0,1) = 
{z E C with \z\ = 1} the series X]n>i is convergent for all t g] 0, 27r[ and by Abel’s 

theorem : 

+ 00 +CX) 

= lim 

r—>-1— 

n=l n=l 

Thus the function m is dehned in £>(0,1)\{1} = {z E C with \z\ < 1, z 1}. By 
positivity we have 

+ 00 +00 

72=1 72=1 

for all 0 < r < 1. Thus for r G [0,1) we get 


+00 

m(r) = Unr"' 

n=l 



ipr 
— ipr 


) 



1 


- dpix) 

— rx 


Lemma 1 For z G T*(0,1)\{1} 


m{z) 


zx 


'[ 0 , 1 ) 


zx 


-dfi{x) 


Proof: First we prove that the function 2 ; 1 -^ -^^^dp{x) is analytic in the open disc 

D{0, 1), continuous in i7(0,1)\{1} and thus is defined. 


Let K denote a compact set in T)(0,This compact is included in £>(0,1)\T)(1,£) 
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with £ > 0 so if G we get zx G -D(0, l)\Zi)(l, e) for all x G [0,1), and thus 

If — zx\ > e. 


Therefore for all z E K and x G]0, 1[ 


zx 

1 — zx 


1 

< -X. 

e 


This proves the continuity of the function jz^dfi^x) over K. Moreover, by 

analyticity ol z ^ ^ ^ [0; 1) the previous boundness condition, the 

function 2 ; 1 -^ j^^dn{x) is analytic in -D(0,1). 

Finally we note that this function coincides with m{r), r G [0,1[, thus by analytic con¬ 
tinuation, we obtain 


m{z) 


zx 


'[ 0 . 1 ) 


zx 


-dfi{x) 


for all z G -D(0,1), and for z = e** with t 7 ^ 0 


+ 00 


mie^^) = 


lim > 


r^l 


Ur,r e 


nAnt 


= lim 


if 

re X 


n=l 


r—>-1— 


'[0,1) 1-reA 


dfiA) = 


if 

ex 


'[ 0 , 1 ) 


1 — e^^x 


dfiA) 


□ 

Extension by continuity of the function -r^ 


Lemma 2 The function can be extended to a continuous function over D{0, 1). 


Proof: The function 1 -|- m doesn’t vanish in -D(0,1) because 


1 -|- m{z) = 


'[ 0 . 1 ) 


1 — zx 


duA) 


and li z = a + ih with -|- 6 ^ < 1 then ax <1 for all x G [0,1) 


so 


Re(l -I- m{z)) = 


1 — ax 


(1 - axY b^x"^ 


dfiA) > 0 


The function is therefore dehned and continuous in -D(0,1)\{1} and analytic in 

D(0,1). 
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In order to study the continuity at the point z = 1, we need to consider two cases : 


1. If is convergent then the series '^UkZ^ is normally convergent in Zi)(0,1) and 


the function m is continuous in -D(0,1) and consequently so is 


l+m ■ 


2. If '^Uk is divergent it is sufficient to prove that \m{z)\ —)> +oo when z ^ 1. Indeed 


m{z) 


1 -^ 1 


l + m{z) IH—^ 

and we extend by 1 at the point z = 1. 

We now show that \m{z)\ —)■ +oo when z —)■ 1. In so doing, consider a sequence of points 
Zk = CLk + ibk with Ofc ^ 1, —)■ 0 and al + bl < 1. Then 


Re(m(2;fc)) = -1 + 


1 — akX 


(1 - ttkxY + blx- 


;d^{x) 


As ttkX < 1 for all integer k and all x G [0,1) the functions 


fk-xe^ 


1 — ttkX 


(1 - afcx)2 + blx'^ 


are positive in [0,1). Using Fatou’s Lemma, we get 


liminf / fkd^{x)> / Imvird fkd^{x) 


'[ 0 , 1 ) 


'[ 0 , 1 ) 


where 


liminf fk{x) = hmfk{x) = - -. 

k 1 — X 


It is therefore sufficient to remark that if ^ Uk is divergent then 


'[ 0 , 1 ) 


1 — X 


dfi{x) > 


> / -- dfi{x) = }^Un = + 00 . 

■'[O,!) ^ ^ n=0 


So for all sequence {zk) converging to 1 in the disc, we have Re{m{zk)) —)■ +cxo and thus 
\m{zk)\ —)■ + 00 . □ 


Convergence of an analytic fnnction T)(0,1), which is continuons on T)(0,1) 
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Lemma [3] provides the use of a Tauberian condition, which turns to be crucial for 
the proof of Theorem 1. 

Lemma 3 Let f{z) = o-n analytic function on -D(0,1), continuous on D{0, 1). 

If f is injective on D{0, 1) then the series is convergent. 

Proof: We proceed in two steps. First, we prove that 



k=l 


On the one hand, the function / is analytic on i5(0,1) thus the image U = f{D{0, 1)) is 
open in C and U is included in the compact set f{D{0, 1)) since / is continuous on £>(0,1). 


On the other hand, the function / being injective on i5(0,1) we have f'{z) ^ 0 for 
all 2 : G L*(0,1) thus / is a diffeomorphism between L)(0,1) and U. By the change of 
variables formula we get 



m 


and with the use of polar coordinates we get the hnite Lebesgue measure of U as the sum 
of the series 


A(f/) = 


n=l 


The convergence of this last series now implies that ^ ^ l^fcl 0. 


Finally, to verify this assertion let > 1 and for n > N write 



k=l 


k=l 


k=N+l 
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Then by Cauchy-Schwarz inequality we have 


1 

n 


k=l 


-l ^ n n 

< _Y^k\h\ + -{ ^ kY'\ ^ k\h\^f/^ 

k=l k=N-\-l k=N-\-l 

N 1 ” 

< -y^k\hu\ + -{ V k\hk\^f/^ 

n ^^ n ^^ 

k=l k=N+l 


and thus 

- n ^ n 

lim sup -V'fc|6fc|<-( ^ k\hk\^Y^‘^. 
Since k -)■ 0 it follows that ^ XlLi k\hk\ 0. 


Taking this intermediate result, we are now in a position to prove Lemma [3l Indeed 
let tn = ^YJk=ikh (with to = 0). We have hn = {K - tn-i) + The series 

X]n>i(^« — tn-i) is convergent because 

N 

— tn-l) = ^ 0, 

n=l 

and thus the series ~ in-i) is Abel summable. Since the series X]n>i 

Abel summable by continuity of /, we get the Abel-summability of the series X]n>i n^n-i- 

Finally, since 

A . 

in—l) ^n—1 t 0, 

n 

it follows from the classical Tauber’s theorem that the series X]n>i n^^-i (and thus 
Zln>i ^n) is convergent. □ 


Proof of Theorem 1 


Using Lemma 2, the function f = a = is dehned and continuous on £>(0,1) and 
analytic in T)(0,1). It follows from Lemma 3 that it remains to prove that / is injective 


in Zi)(0,1). Taking 


m{zi) 

1 + m{zi) 


m{z2) 

1 + m{z2) 


m{zi) 


m{z2), 
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it is sufficient to prove the injectivity of the function m on -D(0,1). Moreover it is sufficient 
to prove this injectivity on -D(0, r) for all 0 < r < 1 : if m is not injective on -D(0,1) then 
there exists Zi ^ Z 2 such that m{zi) = m{z 2 ) ; so m is not injective on D{0,r) where 
r > max(| 2 ;i| , \z 2 \). 

For all 0 < r < 1 the function m is analytic on Zi)(0,r) and continuous on D{0,r). Using 
Darboux’s theorem (Burckel, 1979, p. 310), we could then establish the injectivity of m 
on D{0,r) by showing that m is injective on the circle of radius r with center 0. Indeed, 


let 



We have 



We see that the function (p is decreasing on ]0,7r[ et increasing on jvr, 27r[. It is symmetric 
across tt : we have (p(t) = ip{27r — t). Therefore the only points ti ^ t 2 with (p(ti) = <^(^ 2 ) 
are the pairs {t, 271 — t) with t G [0,7r[. Moreover we have 



hence 


Im(m(re*‘)) = —Im(m(re*^^’^ 


Since we have 



for f G [0, tt), we can’t have 

Im(m(re**)) = Im(m(re**'^"^“*^)) 

Therefore t i-G- m{re^^) is injective. By the Abel Theorem the sum of the series ^ a„ is 
equal to the limit lim^-j.! a(r) where r G [0,1) □ 
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